Observation of the decay <\> — > urn' 
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Abstract 

The reaction e + e~ — > u>ir° 
been studied with SND detector at VEPP-2M 
e + e~ collider in the vicinity of the <fi meson res- 
onance. The observed interference pattern in the 
energy dependence of the cross section is consis- 
tent with existence of the decay 4> — > ujtt with 
a branching ratio of 



B( 



(4.8 



+ 1.9 
-1.7 



±0.8 



10" 



The real and imaginary parts of the decay am- 
plitude were measured. The <ft ^ W7r ° decay was 
observed for the first time. 
PACS: 13.65. +i,14.40.Cs 

1 Introduction 

Recently in experiments with SND and CMD-2 
detectors at VEPP-2M e + e~ collider the study 
of 0-meson rare decays with branching ratios of 
10- 4 -10" 5 became available @, |, |. One of 
such decays is an OZI and G-parity violating 
tfi — > lutt decay. The predicted branching frac- 
tion of this decay is of the order of 5 • 10~ 5 ||. 
It varies within wide limits depending on the na- 
ture of p, uj and </>-meson mixing and existence 
of direct tf> — > lutt transition. Since the nonres- 
onant cross section of the process e + e~ — > unr° 
in the vicinity of the 0-meson peak is relatively 
large, the decay <fi — > uitt° reveals itself as an in- 
terference pattern in the cross section energy de- 
pendence. This allows to extract from the data 
both real and imaginary parts of the decay am- 
plitude. 

The (j> — > u!tt° decay has not been observed 
yet. In our preliminary study |l], ^) an upper 
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limit for the decay probability was placed at a 
level of the theoretical prediction of 5 • 10~ 5 . 

2 Experiment 

The experiment was carried out with SND de- 
tector at VEPP-2M in 1996-1997. SND is a 
general purpose non- magnetic detector Q , based 
on 3.6 t three layer spherical electromagnetic 
calorimeter with 1620 Nal(Tl) crystals cover- 
ing 90% of 47T solid angle. The energy reso- 
lution for electromagnetic showers is ue/E = 
4.2%/ y / E(GeV), the angular resolution is about 
1.5°. The angles of charged particles are mea- 
sured by two cylindrical drift chambers covering 
95% of An solid angle. The angular accuracy of 
charged track measurements is about 0.4° and 
2.0° in azimuth and polar directions respectively. 

Seven successive data taking runs were per- 
formed at 14 energy points in the energy range 
2E=980-1060 MeV. The total integrated lumi- 
nosity AL=3.7 pb' 1 corresponds to 6.3 • 10 6 pro- 
duced (ft mesons. The integrated luminosity was 
measured using e + e~ — > e + e~ and e + e~ — > 77 
reactions with the accuracy of 3%. 

3 Event selection 

For a search for the 4> — > lutt decay we studied 
the cross section of the process 



7r+7r-7r°7r° 



(1) 



in the vicinity of the </>-meson. Events with two 
charged particles and four or more photons were 
selected for analysis. To suppress beam back- 
ground the production point of charged particles 
was required to be within 0.5 cm from the de- 
tector center in the azimuth plane and ±7.5 cm 
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along the beam direction (the longitudinal size 
of the interaction region a z is about 2 cm) . 

Since for the process under study a probability 
to find spurious photons in the events is rather 
large (about 15%), nearly all main <jy meson de- 
cays are a source of background in our search: 



->0 (2) 
K s K L ,K s ^iT + n-- (3) 

-> -> TT+TT^TT ; (4) 



e + e 



e + e — > (f) — > ?7 — > 7T + 7T 7T°. (5) 

The main source of the nonresonant background 
is the process 



7T 7T 7T 7T 



(6) 



which has the same final state as the process (|j) , 
but without the W7r° intermediate state. The in- 
terference between the processes ([[J) and (||) is 
significant only in the small part of the decay 
phase space, where 7r + 7r~7r° invariant mass is 
within the w-meson width. Estimated contribu- 
tion of the interference term in the energy region 
under study does not exceed 2% of the cross sec- 
tion of the process (|l|) and was neglected in the 
further analysis. It was treated instead as an 
additional source of systematic error in the non- 
resonant cross section of the process (|l|). 

To suppress the background from the pro- 
cesses (||) and (||), the following cuts were im- 
posed: 

• average dE/dx losses in the drift chamber 
dE/dx < A-(dE /dx) m i n , where (dE/dx) m in 
are average dE / dx losses of a minimum ion- 
izing particle; 

• the spatial angle between charged particles 
is less than 140°. 

The first condition suppresses events of the pro- 
cess (||) with slow (/3 « 0.25) charged kaons hav- 
ing large dE/dx losses in the SND drift chamber. 
The second condition suppresses events from the 
process (||) with a minimum angle between pions 
from the Ks decay close to 150 degrees. 

For each selected event independent kinematic 
fits were performed under three following as- 
sumptions about the reaction mechanism: 



3. The event is due to the process e + e~ — > 
7r + 7T~7r 7r°. The value of the likelihood 
function Xi-n is calculated together with 
M 37r - the recoil mass of the 7r°-meson clos- 
est to the w-meson mass. 

In the events where the number of photons ex- 
ceeds that required by a certain hypothesis, ex- 
tra photons are considered as spurious ones and 
rejected. To do that all possible subsets with a 
correct number of photons were tested and the 
one, corresponding to a maximum likelihood of 
a certain hypothesis, was selected. 

The distribution of experimental and simu- 
lated events of the process (Eh over the param- 
eter X3n, is shown in Fig. [j]. One can see a 
considerable contribution from the process (|j), 
producing a peak at low X3n- Fig. || shows the 
experimental M rec spectrum, where a clear 77- 
meson peak from the reaction (|5|) is seen. To 
suppress the background from the processes (Q) 
and (||) the following cuts were applied: 

Xz-n > 25, M rec > 620 MeV. 



4 Data analysis 

Figure || shows the distribution of experimental 
and simulated events of the process (EJ) over X4tt ■ 
Considerable difference between the tails of mea- 
sured and simulated spectra indicates the pres- 
ence of a background surviving the cuts. One 
can see that the experimental peak in Fig. ^ is 
broader than a simulated one. It means that the 
simulation of the Xi-n distribution is not precise. 
In Fig.|] the distributions over the 7r° recoil mass 
M^tt for experimental events with Xi-n < 20, sim- 
ulated events of the process @) (the pirn mech- 
anism was assumed), and simulation of the pro- 
cess (fy) with the 11% admixture of the pro- 
cess (|5|) are presented. The last distribution is 
in good agreement with the experimental one. 
For further analysis two additional cuts were ap- 
plied: 

X4* < 40, \M 3n - 782| < 100. 

The total of 4500 events which survived the cuts 
were divided into 4 classes: 



1. 


The event originates from the process 


1. 


X4tv 


< 20, \M 37r 


- 782| 


< 35; 




e + e~ — > 7T + 7r~7r°. The value of the likeli- 














hood function X3n is calculated. 


2. 


X47T 


< 20, |M 3w 


- 782| 


> 35; 


2. 


The event originates from the process 


3. 


X47V 


> 20, |M 3ff 


- 782| 


< 35; 




e + e~ — > 7r + 7r~7r°7. The photon recoil mass 












M rec is calculated. 


4. 


X47T 


> 20, \M 37r 


- 782| 


> 35; 
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with the relative populations of 
0.69:0.09:0.16:0.06 respectively. 

The visible cross section for each class er v i s j 
was fitted according to the following formulae: 



Cvis i 

(E) (E) + /WCE) - 

ow(-E) = £B(u -> 37r)(cr + 
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Da, 



{1 + 5). 



(7) 



Here cr W7r , are the visible cross sections of the 
processes (|l]) and (||) , ctj, ft are the probabilities 
for events of the processes ([!]) or (||) to be found 
in the i-th class (J2 a i — 1 an d ft = 1 at 
B = m^), <70i is the visible cross section of the 
resonant background in the i-th class, <7o is the 
nonresonant cross section of the process e + e~ — > 
urn at E = 77i0, A is its slope, £ is the detection 
efficiency for the process (jl|) at £7 = m^, Z is the 
complex interference amplitude, m^F^, Z?^ = 
— E 2 — iET^(E) are the mass, width and <fi- 
meson propagator function respectively, B(w — > 
37r) = 0.888 is the branching ratio of the ui — > 37r 
decay (5 is a radiative correction calculated 
according to ||. 

The fitting was performed for all 4 classes si- 
multaneously. The class 1 with a small resonant 
background was the most important for the de- 
termination of do, A, Z. Classes 2-4 were used 
to determine the background from <p decays: 
<T 02, o>3, o>4- In the fit it was assumed that for 
the resonant background the distribution over 
is independent of X4tt, and therefore the 
background for the class 1 can be obtained from 
the expression: er^i = (7^2 • (c03/o>4). The va- 
lidity of this relation was checked with the sta- 
tistical accuracy of 15% for simulated events and 
specially selected experimental events of the pro- 
cesses @-(||). The cross section of the resonant 
background o^i in the resonance maximum was 
found to be er^i = (31 ± 17) pb, that is about 4% 
of the visible cross section cr v isi in the class 1. 

The contribution of the process (||) was deter- 
mined using the value of the nonresonant part 
of the cross section in the class 2 and the ra- 
tios aija\ — 0.058 and ft>/ft = 1.05, obtained 
from simulation. The cross section a4 n (E) was 
fitted by a sum of a linear function and an in- 
terference term similar to ^^^(E) in the expres- 
sion (^). For the OZI-rule and G-parity dou- 
ble suppressed <fi — > 7r + 7r~7r 7r decay, the main 
mechanism is a transition of the <f> meson into 
7r + 7T~7r°7r via a virtual photon. The calculated 
real part of the interference amplitude Z due to 
this decay mechanism is equal to 0.127. In the fit 
the real and imaginary parts of the interference 
amplitude of the process (6) were set to 0.127 



and 0.0 respectively, while their errors were set 
to 0.07 and 0.15 in agreement with the accuracy 
of theoretical estimations and experimental data 
on another double suppressed decay <j> — > tt + it~ . 
The nonresonant cross section (||) was found to 
be equal to (11 ± 4)% of the cross section of the 
process ([j]). In the class 1 the contribution of 
this process is smaller: 4.7%. 

Since the simulated distribution over the X4w 
may differ from the experimental one, the co- 
efficients u\ and OL3 were determined from the 
fit. Other coefficients on and ft were derived 
from the relations described above, with the ad- 
ditional assumption that the distributions over 
Xatt in the processes @) and (||) are the same. 
The coefficients o^, ft vary slowly with energy, 
and their energy dependence approximated by 
linear functions was obtained by simulation. 

The detection efficiency e for the process (Q) 
was also obtained from simulation. To estimate 
a systematic error in e, we processed data with 
softer cuts and studied the class of events with 
one charged particle and four photons. It was 
found, that the detection efficiency obtained by 
Monte Carlo simulation must be corrected by 
-11%. The correction is mainly due to the in- 
accuracy in the simulation of distribution 
and loss of charged particles during track recon- 
struction. With this correction, the detection ef- 
ficiency is equal to 17.6 ±1.8% with a systematic 
error of 10%. 

The total number of fit parameters describing 
energy dependence of the cross sections in all 4 
classes of selected events is 11. For each class the 
cross section was measured in 14 energy points. 
At x 2 ' /no = 43/45 the following values of main 
fit parameters were obtained: 

a = (8.2±0.2±0.9)n&, 
A = (0.088 ± 0.009 ± 0.011)nb/MeV, (8) 
Re(Z) = 0.104 ± 0.028 ± 0.006, 
Im{Z) = -0.118 ± 0.030 ± 0.009, 

where the first error is a statistical one obtained 
during the fitting, and the second is systematic 
to be discussed below. 

The visible cross section for the class 1 and 
the fitted curve are shown in Fig. M, Despite 
a 4% resonant background, the interference pat- 
tern is clearly seen. Another representation of 
the interference amplitude is Z — \Z\ • with 



\Z\ = 0.158 ±0.030 ±0.009, 
ip = (-49 ±10 ±4)°. 



(9) 
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The branching ratio for the decay <j> — ► i 
can be obtained from the following relation: 

<Jo ■ \Z\ 2 



B{4> -> coir ) 



0"«) 



(4.8+i;?±0.8)-10" 5 , 
T)/m% = 4220 nb is 



where = 12ttB( 

the cross section in the </>-meson peak 

The following sources of systematic errors 
were considered: inaccuracy of the detection effi- 
ciency estimation, interference between the pro- 
cesses ([j]) and (ft), contribution from the non- 



resonant e + e- 



3n-f reaction, inaccuracy of 



Monte Carlo simulation of the distribution, 
and possible deviations from linear energy de- 



pendence of the nonresonant 



37T7 cross 



section. The total systematic error in the non- 
resonant cross section <jq was estimated to be 
11%. It is determined mainly by the 10% uncer- 
tainty in the detection efficiency. 

Of the factors listed above, the main system- 
atic error in the measured interference amplitude 
Z is introduced by a possible nonlinearity of the 
energy dependence of the cross section of the 
process (|]). It is equal to 4% for Re(Z) and 7% 
for Im(Z). Other factors contribute to Re(Z) 
only at a level of 2%. This nonlinearity also de- 
termines the error in the slope A. 

To estimate systematic errors caused by a 
possible detector instability during lengthy data 
taking, we divided data into three subsets and 
processed them separately. It was found that 
all three subsets are consistent with the fit ob- 
tained from the summary data. The fit param- 
eters obtained from each subset agree well with 
each other and with (||). To check the stability 
of the obtained results, we loosened cuts: for 
the parameter X4tt the cuts 50 and 100 were 
used instead of 20 and 40, for the parameter 
|M 3lr - 782| - 50 and 100 instead of 35 and 100. 
As a result, the detection efficiency increased 
up to 20% and the resonant background in the 
maximum became a^i = (0.20 ± 0.05) nb, that 
is 20% of the visible cross section (||). How- 
ever the obtained interference amplitude was 
Re(Z) = 0.108±0.026, Im{Z) = -0.121 ±0.029 
which does not contradict (|8|), thus confirming 
the validity of the background subtraction pro- 
cedure. As a final result we take (||), because it 
was obtained with a lower resonant background. 
The difference between the results was consid- 
ered as an estimate of the systematic error of 
the resonant background subtraction. Thus, the 
total systematic error of the interference ampli- 
tude Z is equal to 6% for Re(Z) and 8% for 
Im(Z). 



5 Conclusion 



The nonresonant cross section of the process 
e + e~ — > lutt obtained in this work is in agree- 
ment with our old result |lCj on the neutral lo- 
meson decay mode: er(e + e~ — > urn — > 7r 7r°7) = 
(8.7 ± 1.0 ± 0.7) nb. The measured nonres- 
onant cross section ctq is almost twice higher 
than the expected value, calculated with only 
p(770) — > cj7t° transition taken into account. 
The agreement may be significantly improved if 
the known radial excitations of the p meson are 
taken into account. 

The measured interference amplitude Z (^) 
is close to the lower limit of theoretical predic- 
tions although contributions from radial ex- 
citations of p not considered in ||, may affect 
its value. Another important remark is a small 
value of the measured real part of the interfer- 
ence amplitude Re(Z), which could be hardly 
explained by the standard (j> — u> mixing model 
||. For example, the value of 8.2 • 10~ 5 pre- 
dicted in [Q for the branching ratio of the decay 
(j) — > low , is 1.5 times higher than one measured 
in this work. 

The interference amplitude (^) measured in 
this work is five standard deviations above zero. 
Thus, we claim the existence of the decay <p ~> 
u>ir° with the branching ratio of 

B{<f> cjvr ) = (4.8±£? ± 0.8) • 10~ 5 . 
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Figure 2: Distribution over the parameter M rec 
for experimental events and simulation of the 
process e + e~ — > cj7t° — > 7r + 7r~27r° (shaded his- 
togram). The cut is indicated by the arrow. 
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Figure 4: Distribution over parameter for 
experimental events (points with errors) , simula- 
tion of the process (^J) (shaded histogram), and 
simulation of the process (|l|) with an 11% ad- 
mixture of the process (^|) (histogram) . 




Figure 3: Distribution over for experi- 
mental events and simulation of the process ([!]) 
(shaded histogram). 
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Figure 5: Energy dependence of the visible cross 
section for the process (|l]) and optimal curve de- 
scribing data. 



G 



